CHAPITRE 4 



CAS ELEMENTAIRES DE PORTANCE 



Les diffêrents phênomènes de portance en lubrification fluide, peuvent s’êtudier 
sêparêment dans le cas des fluides incompressibles et de viscositê constante sur des exemples 
très simples [1, 2]. Le cas de deux surfaces non parallèles formant un coin d’huile sera êtudiè 
en dêtail au chapitre 5, car il est reprêsentatif de la butêe hydrodynamique. 

1 - CAS DE DEUX SURFACES PARALLELES : PORTANCE HYDROSTATIOUE 



Soit rêcoulement entre deux plaques planes parallèles, de largeur infinie selon Oz. La 
plaque supèrieure de longueur B est fixe ; la plaque infêrieure est animêe d’un mouvement de 
translation uniforme de vitesse = U (fig. 1). Dans ces conditions et compte tenu des 
relations êtablies prêcêdemment, la vitesse du fluide s’êcrit : 
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et l’êquation de Reynolds 



d^ 

dx 2 



= 0 



car la vitesse U et l’êpaisseur h sont 

constantes. Ainsi : — = Cste . 

dx 




Deux cas peuvent se prêsenter : 

a) la pression est identique a l’entrêe et a la sortie. Dans ce cas, dp/dx = 0 ; il n’y a pas de 
portance dans le contact et la rêpartition de vitesse suivant l’èpaisseur du film est linêaire. 
C’est l’êcoulement de Couette. 

b) la pression est diffêrente a l’entrêe et a la sortie. Dans ce cas : 

dp P 2 ~ Pi 

dx L 



ou P[ et P 2 reprêsentent respectivement la pression a l’entrêe et a la sortie du contact. Le 
champ de vitesse dans le film est fonction des pressions imposêes aux extrèmitês du contact. 
Ainsi, la vitesse du fluide est gênêralement due a 2 effets diffêrents : 



l’êcoulement de Couette, c’est le terme 
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l’êcoulement de Poiseuille, c’est le terme 



1 dp 
2p dx 



y(y-h) 
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Dans ce cas il y a portance, mais il faut remarquer que celle-ci est indêpendante de la 
vitesse et de l’êpaisseur du film : c’est une portance hydrostatique. Ces diffêrents types 
d'êcoulements ont êtê prêsentês au chapitre 3. 

2 - CAS DE DEUX SURFACES PARALLELES : EFFET D’ETIREMENT 



Soit l’êcoulement entre deux plaques 
planes parallèles de largeur infinie selon Oz. La 
plaque supêrieure de longueur B est fixe. La 
plaque infêrieure se dêplace en s’êtirant avec 
une vitesse de translation U, = U(x), (fig. 2). 
Ceci est une schêmatisation êlêmentaire du 
formage plasto-hydrodynamique des mêtaux. 

L’êquation de Reynolds s’êcrit : 
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l2 ^ , Fig. 2 : Plaques parallèles : effet d’êtirement 

dp_6qdU 

dx 2 h 2 dx 

soit 

dp _ 6qU 
dx = ~h^~ + 1 
et 

P = J Udx + C[X + C 2 



Les constantes C, et C 2 sont dêterminêes par les conditions aux limites sur la pression. 



Dans le cas particulier ou U = Ax + D et compte tenu des conditions aux limites p = 0 
pour x = 0 et x = B , la pression s’êcrit : 



P 



3qA 

“h 2 ” 
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et la charge pour une largeur L suivant Oz 

W 



Le signe nêgatif de la charge rend compte de l’effet de pompage du fluide qui se 
produit lorsqu’on impose a la distance entre les deux plaques de rester constante et êgale a h. 

3 - CAS DE DEUX SURFACES PARALLELES : EFFET D’ECRASEMENT 



qAB 3 

2h 2 



Soit l’êcoulement entre deux plaques planes parallèles, de largeur infinie suivant Oz. 
La plaque infêrieure est immobile et la plaque supêrieure, de longueur B est animêe d’une 

vitesse V 2 = V (t) suivant Oy , (fig. 3), les deux plaques restent constamment parallèles. 
L’êquation de Reynolds s’êcrit : 
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soit 



et 



d 2 p 12 qV 

3x 2 h 3 

^ = i^x + Cl 

ax h 3 

6 pV 2 

P= r X +C,X + C 9 

h 3 




Fig. 3 : Plaques parallèles : effet d’êcrasement 



Les fonctions C[ (t) et C 2 (t) sont dêterminêes par les conditions aux limites sur la 
pression. 



La rêpartition de pression compte tenu des conditions aux limites p = 0 pour x = 0 et 
pour x = B , s’êcrit : 



P = 



6p V 
h 3 



(x 2 -Bx) 



La portance hydrodynamique W rêsultant de cet effet d’êcrasement s’êcrit : 



W = - 



qVB 3 L 

h 3 



avec 



V = 



dh 

dt 



Dans le cas d’un mouvement d’approche de la plaque supêrieure avec une vitesse de 
module constant |V| suivant Oy, la vitesse est nêgative et l’êpaisseur du film tend vers zêro. 

Pour h > 0 , la portance hydrodynamique due a cet effet d’êcrasement est une fonction du 
temps, elle s’êcrit : 

v|b 3 l 

h 3 




Dans le cas d’un mouvement oscillatoire de la plaque supêrieure autour d’une valeur 
moyenne correspondant a une êpaisseur de film h 0 , l’êquation de l’êpaisseur du film peut se 

mettre sur la forme h(t) = h 0 +hj sin oot ou hj est l’amplitude telle que h t < h 0 et oo la 
pulsation. La vitesse V devient V = hj oo cos oot et la portance hydrodynamique, due a cet 
effet d’êcrasement, s’êcrit : 

LB 3 ph, oo cos oot 

W = 

(h 0 + h[ sin oot) 



a = — avec 0 < a < 1 
h n 



W = 



p oo a B 3 L cos oot 
h 0 (1 + a sin oot) 3 
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soit en posant 




Remarquons que dans ce dernier cas, la pression est positive ou nêgative. L’apparition 
de pressions nêgatives peut entrainer de la cavitation dans le fluide. Par ailleurs il existe un 
dêphasage entre les extrêma d’êpaisseur du film et les extrêma de charge. 

4 - CAS D’UNE DISCONTINUITE 



Soit un patin êchelon, dit patin de Rayleigh (Fig. 4), de largeur infinie selon Oz. II est 
constituê de deux rêgions : une rêgion d’entrêe ou l’êpaisseur du film est constante et êgale a 
h[ , une rêgion de sortie ou l’êpaisseur est êgalement constante et êgale a h 2 (h 2 < hi). Ces 
deux rêgions sont ainsi sêparêes par une discontinuitê dans l’êpaisseur du film, situêe au point 
x = B 0 . Pour êtudier ce problème, il faut isoler la discontinuitê et rêsoudre l’êquation de 
Reynolds sêparêment dans les rêgions d’entrêe et de sortie. 



Dans la rêgion d’entrêe dêfinie par 
O < x < B 0 l’êquation de Reynolds 
s’êcrit : 




soit 

p = Cj x + C 2 




Fig.4 : Patin de Rayleigh 



C, et C 2 sont calculêes a partir des conditions aux limites qui s’êcrivent : 

fp = 0 pour x = 0 

1 p = P m pour x = B 



ou P m est la pression inconnue qui existe au niveau de la discontinuitê. 



II vient ainsi : 



p = — x 

1 r» 



B 



o 



De même dans la zone de sortie dêfinie par B 0 < x < B , l’êquation de Reynolds 



s’êcrit : 



soit 



d 2 p 



= 0 



dx 2 
p = C 3 x + C 4 



Les conditions aux limites 



fP = Pm P our x = B 0 

jp = 0 pour x = B 
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permettent de calculer C 3 et C 4 . II vient : 
P 

(B - x) 



B-B 



o 



Ainsi, la variation de pression est 
linêaire dans le contact (fig. 5). La pression 
inconnue P m est calculêe en êcrivant l’êgalitê 
des dêbits a l’entrêe et a la sortie. 




Fig. 5 - Rêpartition de pression 



Le dêbit pour une largeur L s’êcrit a l’entrêe ou a la sortie : 

0, - L/ 0 h u, dy 



avec 



1 dpi . h ; -y 

u, = L y(y-hi) + — — u 



1 2p, dx L ' h ; 

ou i prend la valeur 1 a l’entrêe et 2 a la sortie. II vient ainsi : 

Lh? dp, LUhj 

t ^ 1 



12 p dx 2 



avec : 



dPi P m 

dx B n 



et 



dp 2 - P m 



dx B - B n 



L’êgalitê des dêbits Q| et Q 2 permet d’êcrire : 



Pm-6^U . h '~ h t 
h l + h 2 



B 0 B-B 0 



soit encore 



P m = 



6pUB s (1 — s) (a — 1) 
h 2 a 3 (l-s) + s 



en posant a = hi/h 2 s = B 0 /B 

pour une largeur L, la charge s’êcrit : 



W = L 



B P„ 



Les valeurs de a et s qui rendent la charge maximale sont 

a = 1,866 ; s = 0,718 
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Fig. 6 : Variation de la charge sans dimension W en fonction du rapport h i/h 2 . 

La figure 6 donne les variations de la charge sans dimension W = h 2 W/ pULB“ en 
fonction de a pour s = 0,718. On peut dêmontrer que pour un encombrement êgal et pour les 
conditions optimales, la charge portante de ce système est la plus êlevêe possible. Cependant, 
les effets d’inertie dans le fluide, lorsqu’ils existent, ont tendance a diminuer la portance du 
patin êchelon. 

5 - APPLICATION 

5.1- Patin a echelon 



On considère un patin a êchelon, ou patin de Rayleigh, semblable a celui schêmatisê 
sur la figure 4. La forme du contact et les conditions de fonctionnement sont telles que : 

Longueur totale : B = 0,5 m 

Longueur de la zone d’entrêe : B 0 = 0,36 m 

Largeur L = 1 m 

Epaisseurs du film : fq =0,2 mm ; h 2 = 0,1 mm 

Vitesse de translation : U = 6m/s 

Viscositê dynamique du lubrifiant : u = 0,01Pa.s 

Si on nêglige l’êcoulement axial, les relations prêsentêes dans le paragraphe 4 
permettent de calculer les caractêristiques de fonctionnement : 

La pression maximale : P m = 1,226 10 6 Pa 
Laportance : W = 306 10 3 N 
Le dêbit : Q = 3,7 10” 4 m 3 /s . 
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La force de frottement est obtenue par intêgration du taux de cisaillement : 



T xy = b 



ÔU 

3 y 



^<2y-h) 



dp 
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-p 
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sur la surface infêrieure : y = 0. On obtient : 



F = -/o L /o B 



(x xy ) y= 0 dx dz = LJ^(h! -h 2 ) + pU 



( B 

h i 



o + B B 0 ^ 



hn 



soit : F = 61,3 + 192 = 253, 3N 

La puissance dissipêe par cisaillement est : 

P = F.U = 1520 W 

et le coefficient de frottement : 

f = — «8,3 10” 4 
W 

5.2 - POMPE A ACTION VISQUEUSE 

On considère le patin a êchelon dêfinit prêcêdemment en 5.1 mais on pratique au 
niveau de l’êchelon une rainure transversale de fa?on a dêriver une partie de Fêcoulement 
(fig. 7). On crêe ainsi une pompe capable de fournir un dêbit Q 0 sous la pression P 0 , si P 0 
dêsigne la pression qui règne au niveau de la discontinuitê. 

Qo , Po 



^IÈIA 




Fig. 7 : Pompe a action visqueuse 



Calculons les caractêristiques (P 0 , Q 0 ) de cette pompe en nêgligeant l’êcoulement 
axial. On distingue 3 rêgions : 

la rêgion d’entrêe I ou l’êpaisseur du film est h , . Le dêbit est donnê par : 



Qi = L 



h 3 p 

n i r o 

12 p B ( 



+ XJ^~ 
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et la force de frottement : 



F, -L 



'Q’o+nuN 

2 hj 



la rêgion centrale III ou la pression est constante et êgale a P 0 . 



la rêgion de sortie II oii l’êpaisseur du film est h 2 . Le dêbit est donnê par 



Qn=L 



hi 



12p, (B - B 0 ) 



■ + U — 



et la force de frottement : 



F„ = L 



-^p 0 + pu^M 

2 0 h. 



La conservation de la masse implique : 

Qi =Qo + Q ii 



d’ou : 



Qo=L - 



12 p 






B 0 B - B 0 



+f[h,-h 2 J 



On remarque que le dêbit Q 0 est une fonction linêaire dêcroissante de la pression P 0 . 
Lorsque le dêbit Q 0 est nul, la pression P 0 est alors identique a la pression P m calculêe 
prêcêdemment en 5.1. Lorsque la pression P 0 tend vers zêro le dêbit Q 0 tend vers la valeur 
limite Q m : 

Q m -Ly<h,-h 2 ) 

Le point de fonctionnement de la pompe (P 0 ,Q 0 ) dêpend des caractêristiques 
hydrauliques du mêcanisme situê en aval de la rainure. Une telle pompe ne fournira jamais un 
dêbit important car la portance de l’êcoulement, conditionnêe par la pression P 0 , serait alors 
très rêduite. 

Considêrons les mêmes donnêes qu’en 5.1 et recherchons les caractêristiques de 

—5 3 

fonctionnement si le dêbit dêrivê Q 0 est de 5-10 ' m / s. (3 f/mn). 



La pression : P 0 « 1,02.10 6 Pa 
la portance : W « 260. 10 3 N 
Lapuissance : P » 1460 W 
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On note que la dêrivation de fluide entraine une perte de portance. Le « rendement » 
de cette pompe est : 

= IjaQo. = 0 035 

P 



C'est un rendement très faible, mais rappelons que l’objectif principal êtait la portance. 
Ce concept de pompe a action visqueuse est quelquefois utilisê pour crêer une autocirculation 
du lubrifiant de la butêe vers le circuit rêfrigêrant [3]. 

Le rendement peut être amêliorê si la portance n’est plus le facteur prêpondêrant. Le 
tableau 1 donne les valeurs adimensionnêes s et a qui conduisent, pour une pression P 0 
donnêe, a un rendement maximal. Ces valeurs sont obtenue pour les donnêes suivantes : 



B h 

B = 0,5m ; L = lm ; U = 6m/s ; h 2 = 10” 4 m ; p = 0,01Pa.s et s = — ; a = — 

B h 2 

Le rendement croit lorsque la pression P 0 diminue, l’êchelon est alors important et la 
rêgion de faible jeu (h 2 ) rêduite. 



Tableau 1 - Evolution du rendement maximal 



P 0 (M.Pa) 


dm 


A 


s 


Q 0 (f/s) 


Qi(f/s) 


Q n (f/s) 


0,3 


0,18 


4,3 


0,94 


0,484 


0,867 


0,384 


0,5 


0,14 


3,2 


0,89 


0,277 


0,653 


0,376 


1 


0,04 


2,1 


0,77 


0,057 


4,29 


0,372 



5.3 - JOINT ANNULAIRE A FUITE 

Un arbre de rayon : R a = 20 mm, animê d’une vitesse de rotation de 3000 tr/mn 
traverse un carter d’êpaisseur : L = 10 mm et d’alêsage de rayon : R c =20,05 mm. A 
Tintêrieur du carter se trouve de l’eau a 20°C et a la pression : P ; = 1 MPa par rapport a la 
pression extêrieure (P e = 0). On suppose que les axes de l’arbre et de l’alêsage sont 
confondus et que Têcoulement est isotherme. Le mêcanisme est schêmatisê sur la figure 7. 

On souhaite connaitre le dêbit de fuite, le couple de frottement et la puissance dissipêe. 

—3 

Pour l’eau, on retiendra une viscositê dynamique : q = 10 Pa.s et une masse volumique : 
p = 10 3 kg/m 3 . 

Le jeu radial C = R c - R a = 0,05 mm êtant très petit devant les rayons 
(C/R a =0,0025) ; on peut nêgliger la courbure du fdm et dêvelopper les surfaces. Le 

problème se ramène a l'êtude de l’êcoulement entre deux surfaces planes parallèles distantes 
de C (fig. 9). Si la surface infêrieure coiTespond a l’arbre, celle-ci est animêe d’un mouvement 
de translation colinêaire a l’axe x, de vitesse U = R a co = 6,28 m/ s . 
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Fig. 8 : Schêmatisation d’un joint annulaire 




L’êcoulement êtant axisymêtrique, les paramètres seront indêpendants de la variable x 
(ou 0 ). Dans ces conditions l’êquation de Reynolds se rêduit a : 




Le champ des vitesses s’êcrit : 



u = -R n (o 

C a 

1 dp , 

w = — — y (y-C) 
2p dz 



et la contrainte de cisaillement x xy a la surface de l’arbre : 
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(’ T 'xy)y=0 M- 



y ) y=0 



= ~h 



Ra M 

c 



Compte tenu des conditions aux limites sur la pression 

jp(z = 0) = P; 

|p(z = L) = 0 



il vient : 



p(z) = Pi| 1-^ 



Le dêbit de fuite axial est donnê par : 

= f 0 f T w(z = L ) Rde d y 



soit : 



Q 



JtRC 3 Pj 
6 p,L 



1,31 10" 4 m 3 /s 



Le couple de frottement est : 

C. - -fofo' R . Vy(y-0)R,d6dr 

soit : 

,, 2jtli R 3 L o) „ , ^ 

C a = — - 3,16 10“ 3 N.m. 

C 

d’ou la puissance dissipêe par cisaillement : 

P = C a oo - 1 W 



Remarques : 

a) La configuration êtudiêe est identique a celle d’un palier lisse centrê. Dans le 
chapitre 7, nous verrons que la position centrêe est une position d’êquilibre souvent instable. 

b) L'effet d'accêlêration du fluide a l'entrêe du joint est nêgligê. Cet effet, qui se 
traduit par une chute brusque de pression a l'entrêe du joint, est appelê effet Lomakin [4]. II 
joue un rôle important sur le comportement dynamique de l’arbre passant a travers le joint. 

c) Les calculs prêcêdents ne sont valables que si le rêgime d’êcoulement est 
laminaire. II faut donc s’assurer que les nombres de Reynolds caractêristiques de l’êcoulement 
sont nettement infêrieurs a 2000 : 
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- nombre de Reynolds de Couette 0f c relatif a la rotation : 



di c = pRa 03 C ^314 
b 

- nombre de Reynolds axial moyen 9i a engendrê par le gradient de pression : 

9t,- pw -C . PQ. ,1043 

p, 2 jt [x R 

Dans les joints la turbulence et souvent due a l'êcoulement axial engendrê par le 
gradient de pression ; cet aspect est dêveloppê dans le chapitre 8. 
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